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We perform N-Body/SPH simulations of disk galaxy formation inside equilibrium spherical and 
triaxial cuspy dark matter halos. We systematically study the disk properties and morphology 
as we increase the numbers of dark matter and gas particles from 10'* to 10^ and change the 
force resolution. The force resolution influences the morphological evolution of the disk quite 
dramatically. Unless the baryon fraction is significantly lower than the universal value, with high 
force resolution a gaseous bar always forms within a billion years after allowing cooling to begin. 
The bar interacts with the disk, transferring angular momentum and increasing its scale length. In 
none of the simulations does the final mass distribution of the baryons obey a single exponential 
profile. Indeed within a few hundred parsecs to a kiloparsec from the center the density rises 
much more steeply than in the rest of the disk, and this is true irrespective of the presence of the 
bar 
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1. Introduction and numerial models 

Although there is a quite extensive literature on the analysis of numerical effects in cosmolog- 
ical SPH simulations of structure formation (e.g. Steinmetz & White 1997), often the complexity 
of such simulations makes it nearly impossible to clearly separate what is physical from what is 
numerical. Attempts have been made to try single out various numerical effects by running sim- 
ple tests on rotating spheres of gas and dark matter and trying to extrapolate results to the case of 
the full cosmological simulations (Okamoto et al. 2003). However, these tests adopt quite ide- 
alised initial conditions that are barely related to the conditions of cosmic structures in a CDM 
model. Moreover, even with such idealised experiments, current resolution tests do not probe a 
large enough range of particle numbers or force resolution to demonstrate convergence. However, 
we wish to understand if we can trust the current state of the art galaxy formation simulations that 
end up of with about 10^ SPH and dark matter particle in single object at z = (e.g. Governato 
et al. 2004). We have studied numerical resolution effects that occur during the dissipative cool- 
ing of gas within isolated equilibrium dark matter halos. At our best resolution we use a million 
particles per halo and hundred parsec force softening. The halo models have structural properties 
and shapes that resemble as closely as possible those arising in CDM simulations. In this high 
resolution regime artificial heating of gas particles by dark matter particles (Steinmetz & White 
1997) should be already under control (A^ > lO'* is required to mitigate those effects). 

In our models we follow the formation of disks inside haloes on scales comparable to the Milky 
Way and M33 (respectively Vdrc ~ 140 and ~ 115 km/s) and with different mass resolutions: LR 
(3 • 10"^), IR (9 • 10"^ ), HR (10^, 5 • 10^) particles in the dark matter and the gas phases respectively. 
Our initial conditions comprise in one case of a spherical CDM like halo (Kazantzidis 2004) built 
with structural parameters expected in the concordance LCDM model for the chosen mass scale, 
and an embedded hot gaseous halo in (approximate) hydrostatic equilibrium. Both components are 
spinning and have the same specific angular momentum distribution, which we take from results 
of cosmological simulations (Bullock et al. 2001). The second set of initial conditions we create 
by merging non-rotating gaseous plus dark matter halos on different orbits (see Moore et al. 2004). 
These produce equilibrium triaxial or axisymmetric prolate and oblate structures with spin param- 
eters in the same range as found within halos from hierarchical simulations. In these paper we 
will show the results from runs using the first type of initial conditions. In our highest resolution 
runs we resolve the gas phase down to IO^Mq or 2 • 10^^% of My,>. The dark matter halos have 
an NFW density disttibution and spin parameter mX = J\E\^ / {GM2 ) from 0.045 to 0. 1 . The gas 
contributes in mass from 6% to 9% of the total mass of the system. We use the parallel Tree+SPH 
code GASOLINE (Wadsley et al (2004)) and a standard cooUng function for a primordial gas of 
helium and (atomic) hydrogen. 

2. Force resolution; disk scalelengths and morphologies 

In this section we will discuss the results of the Milky Way-sized models, which have the 
highest baryon fraction (9%). We ran all simulations with a spUne softening of hi = 0.5 kpc and 
h2 = 2 kpc and focus now on the effects of changing force resolution. The smaller softening leads 
to the formation of a strong gaseous bar of length comparable to the size of the disk at that time, 
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Figure 1: Face on views of the gas density within models having softening hi (left) and hi (right), IR model. 
The upper and lower pictures show a snapshot after 1 Gyrs and 6 Gyrs respectively within a frame of length 
25 kpc on a side. The graph in the right panel shows the logarithmic surface density of the gas disk after 
8 Gyr. Run HR, IR and LR are plotted in red, black and blue respectively. Dashed/solid lines: softening 2 
kpc/0.5 kpc. 



while the disk continues to grow from the inside out with continued gas accretion. Larger values 
of the softening, typical of that used within most cosmological simulations, suppress the bar. Since 
h2 is of order of the disk scale length one expects modes with wavelengths at such scales to be 
damped (with a spline kernel as the one we adopt softening of gravity disappears only at twice /ji). 
In fact in the IR run with softening h2 a bar still forms but dissolves in less than 1 Gyr. 

Assuming conservation of angular momentum and an exponential surface density profile for 
the cold disk material Mo, Mao & White (1998, MMW) calculated the disk scalelength to be 
7?^ = Xr2oo/c °'^/R(^,c,mrf,7£/). The expected scale length for the Milky Way simulations 

should be = 2.2 kpc for X = 0.045 and by replacing r2oo with the cooling radius (which is 
1 10 kpc in HR model) in the simulations. For the HR models we find scalelengths of ~ 2 or ~ 3 
kpc depending on whether the bar is absent or present (smaller scale lengths are found in LR and 
IR models), see Figure 1. The gaseous bar affects heavily the evolution of the baryonic surface 
density, and thus the disk scale length, by transporting angular momentum from the inner part to 
the outer part of the galaxy. This can be seen in Figure ^, where the evolution of the specific angular 
momentum of some disk gas particles which initially formed a ring around the bar, is plotted. Such 
transport of angular momentum of course has a duration comparable to the lifetime of the bar, 
which varies considerably depending on the softening as we noted, and is also affected by the 
adopted mass resolution (number of particles). The transition from a barred galaxy to a nonbarred 
one in the IR model is apparent in Figure 1 and 2. 

An analytic estimate of the stability against bar formation is given by the parameter ( Ef- 
stathiou. Lake & Negroponte 1983, Mayer &Wadsley 2004): = Vpeak I \J {GMdhk /Rd , where 
Vpeak is the peak rotational velocity usually reached around 2/?^ and M^isk is the total mass of the 
disk. We measured for one of our barred simulations at 0.75 Gyr - shortly before it went bar 
unstable - and got £j = 0.60, which is clearly smaller than 0.94,, the value found to be required for 
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Figure 2: Influence of the bar to the evolution of angular momentum of gas disk particles. We followed a 
ring of particles which never reach the inner two kpc, hence they always stay outside the bar: IR model, in 
black: softening 0.5 kpc, in blue: softening 2 kpc, in red: angular momentum of the particles which were 
initially in the bar, softening 2 kpc. How much the angular momentum of such particles grows appears to be 
correlated with the formation of the bar: the bar dissolves for the big softening shortly before 2 Gyr. 




Figure 3: The left panel shows a colour density map of the high resolution M33 gasdisk after 3.5 Gyr. The 
transient spiral patterns closely resemble Sc-Sd galaxies. The right panel shows the surface density profile 
with fitting curves of this model. 



Stability of gas dominated disks against bar formation (see Mayer &Wadsley 2004). 

3. Initial baryon fraction and spin parameter 

The baryon fraction plays an important role. Decreasing the initial baryon fraction to 6% and 
using X = 0.1 lowers the final amount of cold gas by 47% and completely changes the morpholog- 
ical evolution of the disk. Even with high force resolution (softening h = 250pc) the disk remains 
stable and no bar develops, although we found a small nucleus in the centre and some very weak 
spiral arms. 
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This high resolution galaxy model provides an ideal testbed for the analytic model of Van den 
Bosch (2001) for disk structure. Assuming cuspy CDM halos and initial gaseous angular momen- 
tum profiles similar to the dark matter, he finds that it is impossible to form a pure exponential disk, 
in all cases a central luminosity spike results from central low angular momentum material. Indeed 
we are also unable to find a pure exponential disk and we find a steeper profile in the center in our 
simulations (Figure 3). However our M33 model only has a small nuclear gas spike and it is pos- 
sible that with even higher resolution this material collapses to form a dense nucluear star cluster 
as observed within galaxies such as M33. It remains to be seen whether the mass distribution is 
comparable to that in the models by Van den Bosch (2001), which also included star formation but 
do not fully include disk self-gravity. 

We resolve for the first time the mass better than 1O^M0. At such resolution there is clumping 
of infalling gas visible at the edge of the disk; this is the first time that such clumpy gas accretion 
is witnessed in numerical simulations. 

4. Summary 

We simulated the formation of a gaseous disk embedded in a dark matter halo at very high 
resolution. Bar formation is a very common process - particularly with high baryon fractions, low 
spin parameters, and needs high force resolution. Independent of the adopted temperature floor all 
the "Milky Way" disks become bar unstable as soon as we reached a force resolution of 0.5 kpc or 
~ 0.25% of ryir. This picture changes in the M33 model with lower mass and baryon fraction: An 
exponential outer disk plus a small dense exponential nucleus are formed in this case. 

We would Uke to thank SteUos Kazantzidis for providing a code to generate isolated dark 
matter halos. 

References 

[1] Bullock, J., Dekel, A., Kolatt, T., Kravtsov, A., Klypin, A., Porciani, C, Primack, J., 2001, ApJ, 555, 

240 

[2] Efstathiou G., Lake G., Negroponte J., 1982 MNRAS, 199, 1069 

[3] Governato, R, Mayer, L., Wadsley, J., Gardner, J., Willman, Beth, Hayashi, E., Quinn, T., Stadel, J., 
Lake, G., 2004, ApJ, 607, 688 

[4] Kazantzidis S., Magorrian J., Moore B., 2004, ApJ. 601, 37 

[5] Mayer L., Wadsley J., 2004, MNRAS, 347, 277 

[6] Mo, H., Mao S., White S., 1998, MNRAS,295, 319 

[7] Moore B., Kazantzidis S., Diemand D., Stadel J., 2004, MNRAS, 354, 522 

|8] Okamoto T., Jenkins A., Eke, V., Quilis, V., Frenk, C., 2003, MNRAS, 345, 4290 

[9] Steinmetz M., White S.,1997, MNRAS, 288, 545 
[10] Van den Bosch K, 2001, MNRAS, 327, 1334 
[11] Wadsley J., Stadel J., Quinn T., 2004, NewA, 9, 137 



057/5 



